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1. Introduction 

The understanding of gaps-between-jets processes has been subject to great progress 
over the last few years. Of central importance in this context is the energy flow into 
the interjet region as a very useful observable in the description of gaps-between- 
jets processes. According to the method of Sterman et al. [1-4], the cross section for 
interjet energy flow can be 'refactorized' into hard and soft parts at some factorization 
scale /i. For all but the most trivial processes, these parts have a matrix structure in 
the space of colour flows of the hard process, with the hard amplitude represented 
by a vector A in this space and the sum of all possible soft corrections to the cross 
section represented by a matrix S. The \i dependence of each is accounted for by an 
anomalous dimension T, also a matrix in colour space. Logarithms of the ratio of the 
hard and soft scales can be summed to all orders by taking \x of order the hard scale 
in the hard amplitude and the soft scale in the soft matrix, using the exponential of 
the integral of T to connect the two scales. The anomalous dimension matrix T has 
been calculated for all (QCD) 2 — » 2 processes and for various definitions of the final 
state [5-8]. 

It is the purpose of this paper to calculate for the first time the anomalous 
dimension matrix for a 2 — > 3 process, namely qq — > qqg; to distinguish it from the 
one for the 2^2 process qq — > qq, we denote it A. This can serve as a starting 
point to improve the understanding of theoretical aspects in the description of gaps- 
between-jets processes. The calculation of this matrix is also a first step towards 
energy flow analyses of 3-jet processes which are particularly interesting at the LHC. 
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In this paper we consider the process qq — > qqg with a gap defined by a central 
rapidity region of length Y < Ay where Ay is the rapidity separation of the outgoing 
quarks. The real gluon is restricted to the region outside the gap. A is then obtained 
by calculating the virtual corrections to this process from a softer gluon connecting 
the external lines in all possible ways. In common with other calculations of gaps- 
between-jets cross sections [5-8], we assume a perfect real-virtual cancellation outside 
the gap region. Thus the virtual gluon is integrated only over the rapidity interval 
of the gap and over all azimuthal angles. 

We represent the result for A in three colour bases and thereby shed light on 
different aspects of it. 



2. The anomalous dimension matrix for qq — ► qqg 



Let us start with the colour structure. For the qiqj — > quqWa system (where the 
subscripts are the colour indices) there are four independent colour states needed. 
We first choose the t-channel basis 



C 2 = Tm d abc , 

c 4 = T fc ^r bc . 



(2.1) 
(2.2) 
(2.3) 
(2.4) 



The lowest order soft matrix (which contains the traces of the squared operators of 
the basis) is given in this basis by 



S = 



(N c {Nl-l) 


x 










\ 

4) 
N C (N 2 C - 1) 

\N c {Nl -I)) 



(2.5) 



The momenta of the hard process are labeled in the following way 

q(pi) + q(p 2 ) -> q(p 3 ) + <l{Pi) + g(k). 



(2.6) 



We work in the frame in which 1 and 2 collide head on and the gap region is central 
in rapidity, 



Pl = £i(l;0,0, 1), 
p 2 = £ 2 (1; 0,0,-1), 
P3 = Q±3 (cosh y 3 ; 0, l,sinhy 3 ), 
Pi = <?±4 (cosh y 4 ; sin </?, cos if, sinh y±) , 
k — k± (coshy ; sin 0, cos 0, sinh y ). 



(2.7) 
(2.8) 
(2.9) 
(2.10) 
(2.11) 
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Note that in the limit in which the emitted gluon is much softer than the quarks, 
k± <C q±3,4, momentum conservation implies q±3 = q±^ and ip = tt and the kinematics 
are identical to the lowest order process qq — > qq. We are interested in the case that 
the quark jets are either side of the gap and can therefore assume y 3 > and y 4 < 0. 

We denote the rapidity, azimuthal angle and transverse momentum of the virtual 
gluon k! by y', <p' and k' ± , respectively. For future use, we define 

s y = sgn(y). (2.12) 

The gap is defined by a central rapidity region of length Y. Since we are interested 
in the case of k outside the gap and k' within it, we have 

\y'\ <\<\v\ ( 2 - 13 ) 

and hence 

sgn(y - y') = sgn(y) = s y . (2.14) 

We denote the hard amplitude ( [2.6|) evaluated at refactorization scale /i by the 
(four dimensional) vector A(fi). The anomalous dimension matrix A is then defined 
through the evolution of A(fi), 

li4-A= — && (2-15) 
d/x 71 



We can extract A from a one-loop calculation by expanding (|2.15 ) to leading order 



A® = -^ f^AA®, (2-16) 



7T 



where A} ^ and A} 1 ^ are respectively the lowest order and one-loop amplitudes. The 
latter is calculated from the virtual corrections to the hard process from a gluon 
coupling two external lines in all possible ways. We work in the eikonal effective 
theory. The region of integration is 

<</>'< 2tt (2.17) 
-Y/2 < y < Y/2. (2.18) 

Details of the calculation of A can be found in the Appendix. We only state the 
result here written as a sum of four parts. 
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in 

4 ^y 1 




1 v c 



- 271") 






-2s* Y 
4 ^v 1 



- in) - 
—in 



+ 



/ N c \ 
JV C 
iV c 
V N C J 

( C F 








\ 



.(2Y - in) - ^in ) 



\p$\v\) 




(2.19) 



where we have defined 



p(y) = log 



sinh(y/2 + Y/2) 
sinh(j//2 - Y/2) 



Y, 



\ = 



cosh(|y 


+ 


y 


+ Y) - cos(0) 


cosh(|j/| + \y 


-Y)- cos(0) 



-Y 



with 



and 



(p y>0, 
if -(f) y < 0, 



(2.20) 
(2.21) 

(2.22) 



y 



(2.23) 



2/3 y > 0, 
2/4 2/ < 0. 

We have grouped the four terms of A in the following way. For fixed Y, the first line 
contains all the terms that remain when y, y 3 , y A —> ±oo (the high energy limit), the 
second line contains the additional terms that remain when y is finite, the third line 
contains the additional terms that remain when y% or y 4 is finite and the last line 
contains the additional terms that remain when j/3, y 4 and y are finite. For future 
reference we define these four lines to be A 1i2i3i4 respectively. 

For reasons that will become apparent shortly, it will be useful to modify the 
matrix A. Adding a multiple of the identity matrix to any matrix does not change its 
eigenvectors and simply adds a constant to all of its eigenvalues. Moreover, adding 
an imaginary constant to all the eigenvalues of A will not change the physics, since 
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the energy dependence comes from combinations A^ + A^. Therefore we are free 
to add any imaginary constant times the identity matrix to A. From now on we shall 
denote A the matrix obtained from ( [2.19| ) by 



A -»• A + Ay 4 in. 
The eigenvalues of this matrix A are: 

A (1) ~ Y + + N M 2 \V\) + Cf (Jp(2|y s |) + W\y,\)) + ±A, 



2N C 

(2) = Nc y _ Nc + l 
2 2N r 



m + N c \p(2\y\) + C F {\p(2\y z \) + \p(2\y A \)) - ±A, 



A 0) 



N 2 Y - 2m - N C ^N 2 Y 2 - AY in - An 2 

: AN r 



+ N c \p{2\y\) + C F (Jp(2|y 8 |) + Jp(2|y 4 |)) 
A (4) _ N c Y ~ 2i7r + N C ^N 2 Y 2 - AY in - 4tt 2 



N c \ 



(2.24) 

(2.25) 
(2.26) 

(2.27) 



4A"„ 



+ N c \p{2\y\) + C F {\p{2\y,\) + Jp(2|y 4 |)) - f A (2.28) 



and it is diagonalized by 



/ l s _I„ f 3_ 

' 2 v 2 v v 2Af c 



R 



* 2^ + 7T ( 27V C 7T ) S f ^ 



N c 



N c 



N c +2°y N c -2°y 
_ 1 1 

2 2 
1 1 



1 

'2W C 



' 2n 



—s. 



i\< 3 ) 



1 

2JV C 



(2.29) 



/ 



Note that whereas Ai has four different eigenvalues, two eigenvalues of A 4 are degen- 
erate. 

2.1 Block Diagonalization of A 

The anomalous dimension matrix T for the hard process qq — > qq is defined in exact 
analogy to A, (|2.16|) . In the high energy limit (| 2/3,4] — > 00) and in the t-channel 
singlet-octet basis it reads: 







4 4Af| J 



171 



(2.30) 



Two of the eigenvalues of Ai coincide with the eigenvalues of T (to enable this was 
the reason for the modification (|2.24|) of A). We can therefore construct a matrix 



R 



2(N 2 - 1) 



/ 


2 y 


2 v 


Sy 


J-S 

2N C y 


\ 




Nc s 


Nc s 





s y 






N c +2 b y 


N c -2*y 






1 


1 


1 


1 






2 


2 


2N C 




\ 


1 


1 





-1 


/ 



(2.31) 
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which diagonalizes A4 and transforms Ai in the following way: 



R _1 AiR 



fx? 




V 





(2) 












m 



Y 






1 

i_ 



\ 



171 



Z7T 



(2.32) 
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where A]* are the eigenvalues of Aj, which can be obtained from (l2~25H2~28|) by setting 
p and A to 0. Not only is this matrix block diagonal but, remarkably, the upper left 
block is itself diagonal, and the lower right block is identical to T. 
Note that the soft matrix in this basis, 



R f SR 



' 2 



V 



N c +1 
N c +2 









Nl iVc-1 
2 N c -2 

















\ 



(2.33) 



\{Nl-l)J 



is still diagonal and that its lower right block is identical to the soft matrix of the 
2 — > 2 process (the latter property is the reason for our choice of normalization for R). 
It is also interesting to note that in this basis the anomalous dimension matrix A is 
s y independent and that all of the s y dependence is carried by the definitions of the 
basis states, which are different for s y — ±1. 

2.2 The s-channel basis 

Reference [9] advocated using the set of s-channel projectors as the colour basis for 
2^2 processes. In this section we present our results for A in an alternative block- 
diagonal form in which its lower right block is identical to V transformed into the 
s-channel basis and show that its basis states have a simple form. 
For a qq state, the projectors are 



P3 
P 6 



2 (fiki^lj — 8li8kj 



(2.34) 
(2.35) 



We can transform between the t-channel basis we have been using so far and the 
s-channel basis using the matrix 



JVc-l N c +1 



R 



St 



2N C 
-1 



2N C 
+ 1 



(2.36) 



That is, T transforms to 



R st r r. 



N c +1 



si: 



-\{N c -l)Y 



In 



-\{N c + l)Y 
l(N c -l)Y+ N <- 1 



2N C 



I Ti 



(2.37) 
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in the high energy limit and S to 



R S £ S R 



St 



\N C (N C 




\N c {N c +l) 



(2.38) 



Note that the imaginary terms appear only in the diagonal of V and that the entries 
in S correspond to the multiplicities of the basis states, 3 and 6 for N c = 3, two of 
the advantages of the s-channel projector basis. 

We wish to express A in a block diagonal form in which the bottom right block 
is equal to T in the s-channel basis, Eq. (|2.37|) . To this end we define a matrix 



R 



St 



Nc 











\ 



V 










N c 









N c -1 JV e +l 
-1 +1 / 



(2.39) 



in which the bottom right block is equal to R s t and the diagonal entries in the top 
left block are arbitrary - the particular choice made here will lead to a convenient 
result for the soft matrix. 

Transforming A x and S from the original t-channel basis, we obtain 



r/r ^a^rr 



st 



fx? 




\ 





aS 2) 
















\ 



jVc+1 



-\{N c -l)Y 



171 



\(N c + l)Y 



i(jv c -i)y + 



and 



Rg t R^ S R R 



st 



( \N C {N C + 1) 



V 





\N C {N C 








\N C (N C 










N c -1 
2N C 



\ 



171 



(2.40) 



(2.41) 



\N C (N C + 1)J 



By construction, the lower right block of A is equal to T in the high energy limit, the 
lower right block of S is equal to S in the 2 — > 2 s-channel basis and the upper left 
entries of A are left unchanged. The upper left entries of S are set by our arbitrary 
choices in R si for reasons that will be seen shortly. Beyond the high energy limit, the 
bottom right block of A contains a term \C F {p{2\y 3 \) + p(2\y 4 \)) + \N c (p(2\y\) - A) 
times the identity matrix while Y contains just a term ^CF(p(2\y 3 \) + p(2|y 4 |)) times 
the identity matrix, resulting in a small mismatch. 

The actual definitions of the basis states can be read off from the columns of 
RR si and can be written in forms proportional to the s-channel projectors for the 
incoming quarks. Since the matrix R depends on s y , these states are different for 
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Sy — ±1 (recall that s y = +1 implies that the gluon is on the same side of the 
rapidity gap as incoming quark i and outgoing quark k, while s y = — 1 implies that 
it is on the other side). For s y = +1 we have 

Ci oc 
C 2 oc 
C 3 oc 
C 4 oc 

while for s 9 = -1 we 
Ct oc 
C 2 oc 

C 3 oc 
C4 oc 

It is important to note that although we have made arbitrary choices that affect 
the normalizations of these states in order to get S into the form Eq. (|2.41| ), their 
forms are otherwise determined entirely by the physics of A. We see that the two 
states that evolve like a.qq — > qq system have a colour structure given by the qq — > qq 
projectors followed by a gluon emission from the outgoing quark it is on the same 
side of the gap as. The two other states are given similarly by projectors followed by 
an emission from the other outgoing quark, up to colour-suppressed terms coming 
from emission on the same side. 

3. Conclusions 

We have calculated the anomalous dimension matrix A for the five-parton process 
qq — > qqg and presented it in several different colour bases. It seems likely that the 
generalization of the s-channel basis, Eq. ( |2.4U| ), will be most useful both for gaining 
insight into the physics of A and for performing practical calculations. We anticipate 
using A to improve the theoretical understanding of gaps-between-jets processes and 
ultimately to calculate energy flow observables in 3-jet processes, which are particu- 
larly interesting at the LHC The latter however requires the anomalous dimension 
matrices for all 3-jet processes to be calculated, a highly non-trivial problem: in the 
most complicated case of gg — > ggg one expects up to 44 independent colour ampli- 
tudes and a deeper theoretical insight seems necessary to organize the calculation. 



2 (^mi^nj ^ni^mj 



have 



2 (^rni^nj $ni$mj 

^{s mi s nj - s ni s mj 

2 (^mi^nj ^ni^mj 



bkmTln _|_ ^km^ln ] ■ 

bkmTln + jy _ ^krrfiln ) ■ 



1 



N c + 1 



-bkmXln j , 



N c -1 



-bkmTi n I , 



^krnXln j , 



SkmT^ ) . 



(2.42) 

(2.43) 
(2.44) 
(2.45) 

(2.46) 
(2.47) 
(2.48) 
(2.49) 
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In particular it would be extremely interesting to see whether the block diagonal 
structure we found for qq — > qqg can be generalized to arbitrary processes, with one 
block always equal to the anomalous dimension matrix of a lower-order process. 



A. Calculation of A 



In the basis 



A has the following colour structure 



/ 5^(012 + n 34 + n 14 + n 23 ) (i - ^)(n 12 + n 34 + n u + n 23 ) 



n 12 + n M + n 14 + n 23 -^(n^ + + n 14 + n 23 ) 



A 



+f (fi 14 + n 23 ) 
o 



(f -j_)(-n 14 + n 23 ) 



^-(Hl2 + ^34 + 0l4 + ^23) 



i(-f2 14 + n 2:i ) 
f (-n M + n 23 ) 

K-Qis + Om) 



-n 12 + n 34 — 2^-(n 12 + n 34 + n 14 + n 23 ) 



V 



+ 



+ 



+f (n M + a») 



/ 



o 





o f (n u . + n 2 , + n 3 , + n 4fc ) 

\ n lk + n 2k + n sk + n ik (^-±)(n lh + n 2k + n 3k + n ik ) -n lk + a 2k - n 3k + f(n lk - u 2k - n 3k + n ih ) J 

n ^ri-<» l; ■ !>,:! \ 



^ ^(f^lft — ^2* — ^3* + f^4fc) 


^■(S^l/fc + — ^3* — ^4k) 



/(f -2^)^13 + ^24) 



f (-n 13 + n 24 ) 
\ o 



f (-n 13 + n 24 ) o 
-^(n 13 + n 24 ) o o 

(f " 2^)^13 + M o 

"2^(^13 + ^4)7 



(A.l) 



Qij corresponds to the case in which the virtual gluon couples quarks i and j. 
accounts for the coupling of quark i and the gluon k. These functions are given by 



and 



n 



ik 



-SiS k Sg 



[ dy'ir- u 'ik - 5(1 - <W i7r 
.J si 2n 



where we have introduced the shorthands 



\k'1 pj-pj 
Pi-k' k'-pj 



ik 



2 Pi ' k 

Pi-k' k'-k 



(A.2) 



(A.3) 



(A.4) 



We have 6^5 j = — 1 if i and j are both incoming or both outgoing and +1 otherwise. 
S g depends on the topology of the triple-gluon vertex: in our convention in which the 
indices of if abc are labeled in an anticlockwise direction around the vertex, if, with 
the vertex rotated so that the momentum of the eikonal gluon is flowing horizontally 
from left to right, the soft gluon is above it, then 5 g = +1, and if below, 5 g = —1. 

It is worth pointing out that in the general case, the elements of A, Eq. (|A.1| ) 
obey the following equality: 



= ( n ° SUm ) 



(A.5) 
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implying that in an orthonormal basis in which S is equal to the identity matrix, A is 
a symmetric matrix. This property is therefore valid independently of the observable 
to which A contributes. It was pointed out in Ref. [10] that this property is true 
of all anomalous dimension matrices that have been calculated to date, although no 
explanation of this fact was offered. 

Carrying out the integrations in ( |A.3| , |A.4| ) over the region Q: 

<</>'< 2vr (A.6) 
-Y/2 <y'< Y/2, (A.7) 



we obtain 



£ll2 
^34 

^23 
^3fc 



-l(Y-m)-\p(2\y 3 \) 
l -Y + l -p{2\y,\\ 
ly + ip(2|y 3 |), 



1 



P(2|2/4|), 



p(2|y 3 |), 



P(2|2/ 4 |), 



1 



s y )Y + -p(2\y\), 



-\(l + s y )Y-±p(2\y\), 



n 



Ik 



where we have defined 



l(l-s y )Y + ~p(2\y 3 \) + lp(2\y\ 



--(l + Sy )\(\y 3 \ + \y\ 



Z7T 



\^ + s y )Y+ l - P {2\y,\) + l -p{2\y\) 

-^i 1 - s v)H\yi\ + \vl<p - f 



ITi 



KvA) 

It is useful to note that 



s y = sgn(y), 

snMy/2 ±m _ 
P[V) ~ g sinh(y/2 - F/2) ' 

1 cosh(y + Y) — cos(0) 
— log 

2 cosh(?/ — K) — cos(0) 



— Y. 



(A.8) 
(A.9) 
(A.10) 
(A.ll) 
(A.12) 
(A.13) 
(A. 14) 
(A.15) 

(A.16) 



\(y,0) = p(y). 



(A.17) 

(A.18) 
(A.19) 

(A.20) 
(A.21) 



- 10 - 



Acknowledgments 

We are extremely grateful to Jeff Forshaw for interesting discussions of this and 
related topics. AK is supported by PPARC grant PP/D000157/1. 

References 

[1] J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 308, 833 (1988). 

[2] M. G. Sotiropoulos and G. Sterman, Nucl. Phys. B 419, 59 (1994) 

[3] H. Contopanagos, E. Laenen and G. Sterman, Nucl. Phys. B 484, 303 (1997) 

[4] N. Kidonakis, G. Oderda and G. Sterman, Nucl. Phys. B 525 (1998) 299 

[5] G. Oderda and G. Sterman, Phys. Rev. Lett. 81, 3591 (1998) 

[6] G. Oderda, Phys. Rev. D 61 (2000) 014004 

[7] C. F. Berger, T. Kucs and G. Sterman, Phys. Rev. D 65 (2002) 094031 
[8] R. B. Appleby and M. H. Seymour, JHEP 0309 (2003) 056 
[9] Yu. L. Dokshitzer and G. Marchesini, arXiv:hep-ph/0508130. 
[10] M. H. Seymour, arXiv:hep-ph/0508305. 



- 11 - 



